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CYCLES 



ASHKAN NIKEGHBALI AND DIRK ZEINDLER 



Abstract. The goal of this paper is to analyse the asymptotic behavior of the cycle process 
and the total number of cycles of weighted and generalized weighted random permutations 
which are relevant models in physics and which extend the Ewens measure. We combine tools 
from combinatorics and complex analysis (e.g. singularity analysis of generating functions) 
to prove that under some analytic conditions (on relevant generating functions) the cycle 
process converges to a vector of independent Poisson variables and to establish a central 
limit theorem for the total number of cycles. Our methods allow us to obtain an asymptotic 
estimate of the characteristic functions of the different random vectors of interest together 
with an error estimate, thus having a control on the speed of convergence. In fact we 
are able to prove a finer convergence for the total number of cycles, namely mod-Poisson 
convergence. From there we apply previous results on mod-Poisson convergence to obtain 
Poisson approximation for the total number of cycles as well as large deviations estimates. 
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1. Introduction 



In this paper we are interested in finding the asymptotic behaviour of the cycle structure 
and the number of cycles of weighted random permutations which appear in mathematical 
biology and in theoretical physics. More precisely, we define the weighted and generalized 
weighted probability measures on the group of permutations Sn of order n as follows: 

Definition 1.1. Let a & Sn be given. We write Cj{a) for the number of cycles of length j 
in the decomposition of cr as a product of disjoint cycles (see also Definition 2.2). 

(1) Let 9 = {Om)m>i given, with 9j > for every j > I. We define the generalized 
weighted measures as 



H := ^ n (1-1^ 



m=l 



with hn = hn{Q) a normalization constant and Hq := 1. 
(2) More generally, we define the generalized weighted measures as follows. Let '■ 
N — )■ M>o be given for m > 1 with Fm{0) = 1. We then define 

1 " 

^^[^]^=;i!M^n^™M^)) (1-2) 



m=l 



with hn{F) a normalization constant. It follows immediately from the definition that 
Pe[.] =Pf[.] with F^{k) = el. 

When 9j = 1 for all j > 1, the measure Pe [•] is the uniform measure on permutations and 
these are well studied objects with a long history (see e.g. the monograph [1] for a complete 
account). In particular it was proven by Goncharov ([IZ|) and Shepp and Loyd ([26]) that 
the process of cycle counts converges in distribution to a Poisson process, that is as n — )■ oo, 

C[''\ct\..) ^'iZ,,Z,,...), (1.3) 

where the Zj are independent Poisson distributed random variables with ^[Zj] = -. There 
have also been further studies to analyze the rate of convergence for 



P 



where 6 is a fixed integer. The above can be interpreted as a result on small cycles. There 
exist as well results on large cycles, due to Kingman ([22]) and Vershik and Shmidt (j27j) 
who prove that the vector of renormalized and ordered length of the cycles converges in law 
to a Poisson-Dirichlet distribution with parameter 1. Moreover if one notes 

Ron = Cf") + . . . + C(") 

the number of cycles, then the distribution of is well known and a central limit theorem 
can be proven: 



- logn 
2 



^'Ar(0,l), (1.4) 



where A/'(0, 1) stands for a standard Gaussian random variable. In fact one can prove Poisson 
approximation results for i^on — 1 as well as large deviations estimates. For instance Hwang 
(|19j) showed that for k ~ xlogn, 



TOFT^ u (log^)^ ^exp(-logn; , . 
P[Kon = k\ = — — — — - + O 



{k-l)\ Vr(l + r) V(logn^ 

where r = {k — l)/logn. 

Similar results exist if one considers the more general Ewens measure corresponding to 



9j = 9 > for all j in equation (1.1) defining Pe [.]. This measure was introduced by Ewens 
|10 ) to study population dynamics and has received much attention since. In particular 



(1.3) and (1.4) hold with ^[Zj] = j and logn replaced with 6'logn in the central limit 



theorem. Estimates on the rate of convergence as well as Poisson approximation results for 
Kon area available as well. 

The measure Pe [.] is thus a natural extension of the Ewens measure and besides has a 
physical interpretation: indeed such a model appears in the study of the quantum Bose 
gas in statistical mechanics (see [7], [6] and [5]). There it is of interest to understand the 
structure of the cycles when the asymptotic behaviour of 6j is fixed. The case where 6j — )■ 6, 
i.e. asymptotically the Ewens case, was also considered in [3]. Another random variable of 
interest that we shall not consider in this paper is Li, the length of the cycle that contains 
1, which can be interpreted as giving the length of a typical cycle, has also been considered 
in some of the above mentioned papers. It appeared in these works that obtaining the 
convergence in distribution of the cycle process or the central limit theorem for the number 
of cycles is a challenging problem. Indeed, there does not exist something such as the 
Feller coupling for the random permutations under the measure Pe [.] , since in general these 
measures do not possess any compatibility property between the different dimensions. The 
main important property of Pe [•] is that it is invariant on conjugacy classes, and we shall 
exploit this fact. 

In a recent skillful paper, Ercolani and Ueltschi ([H]) have obtained, under a variety of 
assumptions on the asymptotic properties of the the convergence of the cycle process 



to a Poisson vector, as in (1.3), with this time ^[Zj] = -j-. In some cases they obtained an 
equivalent for E[i^'on] and in some "degenerate cases" they proved that the total number of 
cycles converges in distribution to 1 or to 1 plus a Poisson random variable. But their method 
which is a subtle saddle point analysis of generating functions does not give any information 
on the asymptotic behaviour of the different random variables under consideration nor on 
the rate of convergence. But on the other hand because the method is general, they are able 
to cover concrete cases corresponding to a large variety of assumptions on the asymptotic 
behaviour of the 6'j's. 

The goal of this paper is to bring a complementary point of view to the approach of Ercolani 
and Ueltschi ([9]) and to provide sufficient conditions on the 6'j's, or more precisely on the 



analyticity properties of the generating series of {9m)m>i-, under which one has (1.3), (1.4), 
estimates on the rate of convergence, as well as Poisson approximation and large deviations 
estimates for the total number of cycles K^n- 

Our approach is based on the so called singularity analysis of the generating series of (6'„)„>i 
and is general enough to deal with the case of the more general measures Pf [.]. The starting 
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point is the well known relation 

oo oo Q 

J] h^e = exp {g{t)) , with g{t) = f ^" 

n=l n=l 

which relates the generating series of the sequence (/i„)„>i to that of the sequence (6'„)„>i, 
and its generalized version for the measure Pp [.] (in fact these formulas will follow from 
some combinatorial lemmas which are useful to compute expectations or statistics-e.g. the 
characteristic function- under Pp [.]). To obtain an asymptotic for /i„ as well as an estimate 
for the characteristic functions of the different random variables under consideration, it will 
reveal crucial to extract precise asymptotic information with an error term on the coefficient 
of w)] for F(t,w) = exp{wg(t))S(t,w) where S{t,w) is some holomorphic function 

in a domain containing {{t,w) G C^; |t| < r, |iy| < f} where r is the radius of convergence of 
the generating series g{t). This will reveal possible if one makes some assumptions on the 
analyticity property of the generating series g{t) together with assumptions on the nature 
of its singularities at the point r on the circle of convergence. Several results of this nature 
exist in the literature (see e.g. the monograph [T^) but the relevant ones for us are due to 
Hwang ([20] and [18j). Combining these results with some combinatorial lemmas, we are 
abl e to show that the cycle count process converges in law to a vector of Poisson process as 



in (1.3), with this time ^[Zj] = where we recall that r is the radius of convergence of 
the generating series g(t) of the sequence (6'„)„>o. We also have a an estimate for the speed 
of convergence. We are also able to prove a central limit theorem for K^n, as well as Poisson 
approximation results and large deviations estimates. In fact our methods allow us to prove a 
stronger convergence result than the central limit theorem, namely mod-Poisson convergence. 
This type of convergence together with mod-Gaussian convergence was introduced in [21j 
and further developed in [23]. It can be viewed as a higher order central limit theorem from 
which one can deduce many relevant information (in particular the central limit theorem). 
More precisely the paper is organized as follows: 

• In Section [2] we establish some basic combinatorial lemmas and review some facts 
from complex analysis (e.g. singularity analysis theorems and the Lindelof integral 
representation theorem) from which we deduce the asymptotic behaviour of the nor- 
malization constant hn (under some assumptions on the generating function g{t)). 

• In Section [3] we prove the Poisson convergence for the cycle process together with a 
rate of convergence; 

• Section 4 is devoted to various limit theorems for the total number of cycles K^n'-, 

• Section 5 contains some examples; 

• In Section [6] we prove general limit theorems under the more general measure Pp [.] . 
We shall illustrate these results with the example of exp-polynomial weights and a 
toy example of spatial random permutations which plays an important role in physics 
(see e.g. [3]). Here we consider the simpler case where the lattice is fixed. It is our 
hope that our methods can be adapted to deal with more complicated cases where 
the lattice is not fixed anymore. 

2. Combinatorics and singularity analysis 

2.1. Combinatorics of Sn and generating functions. We recall in this section some 
basic facts about Sn and partitions, and at the end of the section state a useful lemma to 



perform averages on the symmetric group. We only give here a very short overview and refer 
to [T] and [21] for more details. 

We first analyse the conjugation classes of Sn since all probability measures and functions 
considered in this paper are invariant under conjugation. It is well known that the conjuga- 
tion classes of Sn can be parametrized with partitions of n. 

Definition 2.1. A partition X is a sequence of non-negative integers Ai > A2 > ■ ■ ■ eventu- 
ally trailing to O's, which we usually omit. We use the notation A = (Ai, A2, ■ ■ ■ , A;). 
The length /(A) of \ is the largest I such that \i 7^ 0. We define the size |A| := Ylim '^m- We 
call A a partition of n if |A| = n. We use the notation 

J2i-)--= E {■■) and ■■= E (••)• 

AHn A partition of n A A partition 

Let 0" G S'„ be arbitrary. We can write a = ai - ■ ■ ai with (Xj disjoint cycles of length Aj. 
Since disjoint cycles commute, we can assume that Ai > A2 > ■ ■ ■ > A^. We call the partition 
A = (Ai, A2, ■ ■ ■ , A;) the cycle-type of a. We write C\ for the set of all cr G 5"^ with cycle 
type A. One can show that two elements cr,T & Sn are conjugate if and only if a and r have 
the same cycle-type and that the Cx are the conjugacy classes of Sn (see e.g. [21] for more 
details) . 

Definition 2.2. Let a E Sn be given with cycle-type A. The cycle numbers Cm and the total 
number of cycles K^n are defined as 

n 

Cm = := # A. = m} and Kon ■= ^ ^^^^ ■ (2-1) 

m=l 

The functions C^ and K^n depend only on the cycle type and are thus class functions (i.e. 
they are constant on conjugacy classes). 

All expectations in this paper have the form ^ Xlo-es '^i'^) for a certain class function u. 
Since u is constant on conjugacy classes, it is more natural to sum over all conjugacy classes. 
We thus need to know the size of each conjugacy class. 

Lemma 2.3. We have 

\S I " 

\Cx\ = with zx ■■= TT m'^'^Cm! and Cm = Cm(A) = # {A,; A^ = m} , (2.2) 

i ! 1 I 

and 

U&Sn A 

for a class function u : Sn ^ C 

Proof. The first part can be found in [21] or in [S] chapter 39]. The second part follows 
immediately from the first part. □ 

Given a sequence (5'n)n>i of numbers, one can encode the information about this sequence 
into a formal power series called the generating series. 
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(2.3) 



Definition 2.4. Let {gn)„^f,j be a sequence of complex numbers. We then define the (ordi- 
nary) generating function of {gn)^^^^ (is the formal power series 

oo 

Git)=G{gn,t) = J29nt''- (2-4) 

n=0 

Definition 2.5. Let G{t) = ^^o^'"'^" formal power series. We then define [t"] [G] : = 
gn, i.e. the coefficient oft"- in G(t) 

The reason why generating functions are useful is that it is often possible to write down 
a generating function without knowing Qn explicitly. Then one can try to use tools from 
analysis to extract information about gn, for large n, from the generating function. It should 
be noted that there are several variants in the definition of generating series and we shall 
use several of them and still call all of them generating series without risk of confusion. We 
will also later replace (see Section |4]) gn by holomorphic functions gn{w)- Such generating 
functions are then called bivariate generating functions. Again for simplicity we shall still 
call them generating functions. 

We now introduce two generating functions (in the broad sense) which will play a crucial role 
in our study of random weighted permutations under the measure Pe [(•)]• Fo^' © = (^n)n>i, 
we set 

oo „ / °^ \ 

■= E Geit) := exp E (2.5) 

fc=i \fc=i / 



For now, gQ{t) and GQ{t) are just formal power series. We will see in Section |3] and in 
Section ^ that the asymptotic behaviour of Gm^ and Kon depend on the analytic properties 
of gelt), essentially because of the remarkable well known identity (that we will quickly 
derive below) 



Ge(t) = X]M". (2.6) 



n=0 



One of the main tools in this paper to compute generating series is the following lemma (or 
cycle index theorem) of which we shall prove a more general form in Section |6] to deal with 
the more general measure Pp [.] . 

Lemma 2.6. Let (am)meN be a sequence of complex numbers. Then 

= E i (n («.»«"■)"•") = -P (e l^a^r) (2.7) 

A ^ \m=l j A ^ \m=l / \m=l / 

with the same z\ as in Lemma \2. 3[ 

If one of the sums in (2.7) is absolutely convergent then so are the others. 

Proof. The first equality follows immediately from the definition of Gm- The proof of the 
second equality in (2.7) can be found in [21] or can be directly verified using the definitions of 
zx and the exponential function. The last statement follows from dominated convergence. □ 
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We now use this lemma to show (2.6), i.e. G@{t) = Yl^=o^'n.t^- know from (1.1) that 

creSn m=l 

It now follows from lemma 12.61 that 



(2.9) 



Ahn ^ m=l 



Ahn ^'^ m=l 



oo ^ IW / CO „ \ 

ra=l A m=l \m=l / 



This proves (2.6 ). 



2.2. Singularity analysis. If a generating function g{t) is given then a natural question is: 
what is [t"] [g\ and what is the asymptotic behaviour of [t"] [g\. If (/(t) is holomorphic near 
then one can use Cauchy's integral formula to do this. Unfortunately it is often difficult to 
compute the integral explicitly, but there exist several other results which allow to achieve 
this task. One such theorem is due to Hwang [18], and we prepare it with some preliminary 
definition and notation. 



Definition 2.7. Let < r < R and < < | be given. We then define 

Aq = Ao{r,R,(f)) = {z eC; \z\ < R,z ^r,\ aig{z - r)| > 0} 



(2.10) 



Figure 1. Illustration of Aq 

Definition 2.8. Let g(t) and 6 > 0,r > be given. We then call g{t) of class J^{r, 6) if 

(1) there exists R > r and < < | such that g{t) is holomorphic in Ao(r, R, 0), 

(2) there exists a constant K such that 

git) =e\og{ — ^—— \ +K + 0{t-r) for t ^ r. (2.11) 



1-t/r 

We shall use the shorter notation g(t) G J^{r, 



We emphasize only the dependence on 9 and r since the other constants do not appear in 
the main results. We now state an important theorem due to Hwang: 

Theorem 2.9 (Hwang [ig). Let F{t,w) = e'^sd) s{t,w) he given. Suppose that 

(1) Let g{t) is of class J^{r, 6), 

(2) S(t,w) is holomorphic in {t,w) for \t\ < r and \w\ < r for some r > 0, i.e. S{t,w) is 
holomorphic in {t,w) in a domain containing the set {{t,w) e C^; \t\ < r, \w\ < f}. 

Then 



uniformly for \w\ < r and with the same K as in (2.11) 



Proof. The idea of the proof is to take a suitable Hankel contour and to estimate the integral 
over each piece. The details can be found in [HI chapter 5]. □ 

Remark: 

We use most times w = 1 and S independent of w. One can also compute lower order error 
terms if one has more terms in the expansion of g{t) near r. 

A natural question at this point is: how can one prove that g{t) is of class -F(r, 9)1 It is most 
times easy to compute the radius of convergence of g{t), but it is not obvious how to show 
that g{t) is holomorphic in some Aq. A way to achieve this is through Lindelof's integral 
represention: 

Theorem 2.10 (Lindelof's integral represention). Let (f){z) be a holomorphic function for 
Re{z) > 0, satisfying 

|0(z)| < Ce^l^l for \z\ oo andRe{z) > ^ with some A G]0,7r[,C > 0. (2.13) 
Let g(t) := YlT=i 4>{k){—t)'' . The radius of convergence of g(t) is e""^ and 

g{t) = — / 0(^)t^_^ dz. (2.14) 

Furthermore g{t) can be holomorphically continued to the sector —{ti — A) < arg(t) < (tt — A). 
Proof See [121 Theorem 2]. □ 



In many situations Theorem 2.10 allows us to prove holomorphicity in a domain Aq, but 
does not give any information about the behaviour of g{t) near the singularity. One way to 
compute the asymptotic behaviour of g{t) near the singularity is to use the Mellin transform. 
We do not introduce here the Mellin transform since this would take us to far away from the 
topic of this paper. We would rather refer to [13] for an introduction and to [HI Section 3] 
for an application to the polylogarithm. 



Theorem 2.10 is not always so easy to apply and the computation of the asymptotic behaviour 
near the singularity is often very difficult. An alternative approach is to combine singularity 
analysis with more elementary methods. The idea is to write F = F1F2 in a. way that we 
can apply singularity analysis to Fi and can estimate the growth rate of [t"] [F2]. One can 
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then compute the coefficient [t"] [F] directly and conclude with elementary analysis. This 
method is called the convolution method. We now introduce a class of functions which is 
suitable for this approach. 

Definition 2.11. Let g{t) and 6'>0, r>0,0<7<l he given. We then call git) of class 
eJ^{r, 9, 7) if there exists an analytic function go in the disc {\z\ < r} such that 



g{t) = eiog 



+ go{t) fort^r with [r] [^0] = 0(r-"n-^-^) 



1-t/r 

We can now state a singularity analysis theorem corresponding to the class eJ-'(r, ^,7): 

Theorem 2.12 (Hwang [20]). Let F{t,w) = 6""^^^^ S(t,w) he given such that 

(1) g{t) is of class eJ-'(r, 6*, 7), 

(2) ^(t, w) is holomorphic in {t, w) for \t\ < r and \w\ < r for some f > 0. 
We then have 

S{r, w) 



(2.15) 



r] [F{t,w)] 



+ Rn{w) 



(2.16) 



with K = go{r) and 



RJw) 



O ( ^ 



,^9Rc(,j,)-l-7jQg(„) 



, ifRe{w)>0 
ifRe{w) < 0. 



uniformly for hounded w. 
Proof See [20]. 

We can now compute the asymptotic behaviour of 

Lemma 2.13. Let ge{t) he of class J^{r,9) or of class eJ^{r,9,-y). We then have 

1 m 



(2.17) 



□ 



— , n — )■ 00. 



(2.18) 



with the same K as in ( 2.11^ resp. as in Theorem 2.12 



Proof. We have proven that J2'^=o^nt^ = Ge(t) = exp{ge{t)). We thus can apply Theo- 
rem 



2.9 resp. Theorem 2.12 with g(t) = ge{t) for w = 1 and S{t, w) = 1. We get 

hr 



r"r(0) 



n — )■ 00 



(2.19) 

□ 



Remark: 

There exist several other versions of the Theorems 2.9| and Theorems 2.12 For example one 
can replace log(l — t/r) by other functions (see [Hj) or allow more than one singularity (see 
[151 chapter VI. 5]. 
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3. Limit theorem for the cycle numbers 



In this section we establish the convergence in distribution of the cycle process to a vector 
of Poisson random variables. 

Theorem 3.1. Let b be fixed. We then have as formal power series 



n=0 



m=l 



If go is of class J-'{9,r), then 

b 



exp I 



n) 
m 



m=l 



exp 



vm=l 



\m=l 



m 



o 



(3.2) 



If ge{t) is of class eJ^{r,9,'y), then 
Ee 



exp ( z ^ s^C^A = exp [ ^ ^(e^™ - l)rA +o( 

\ m=l / J \m=l / ^ 



(3.3) 



The convergence result now follows immediately from Theorem 3.1 



Corollary 3.1.1. Let <d = {dmjmen be given and Sn be endowed with Fq [(.)]. Assume that 
ge{t) is of class Ti^r, 6) or of class eJ^{r, 6, 7). We then have for each 6 G N 



(3.4) 



with Yi, - ■ ■ ,Yb independent Poisson distributed random variables with E [1^ 
Let Ci, ■ ■ ■ , c„ be integers. We then have for g^ G J^(r, 6) 



Cf"^ =ci,--- ,C,("^ = CfcJ -Pe [Fi = ci,--- ,n = Cb] 
and for g^ G eJ-'(r, 6', 7) 

fcf") = ci, ■ ■ ■ , C^") = cj - Pe [Fi = ci, ■ ■ ■ , n = c^] 



O 



log(n) 



Proof. The first part follows immediately from Levy's continuity theorem. To prove the 
second part, we use the Fourier inversion formula. Let ?/'(si,-- - , sj,) be the characteris- 
tic function of c{"\ C^"\ • • • and 0(si, ■ ■ ■ , Sb) the characteristic function of Fi, ■ ■ ■ ,1^6. 

with 0(.) 



Equation (3.2) resp. (3.3) shows that — (t)\ = O (^) resp. — = O 
independent of Si, ■ ■ ■ , s^. On the other hand we have 



log(") 



P« 



1 



(27r)^ 

This proves the corollary 



(V'(si,-- - ,S6) -0(si,--- ,Sfe))e"^™=i^"™''" dsi---dsb. 



□ 
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Proof of Theorem 3J_. We first compute the generating function in (3.1). Tlie factor hn in 
(3.1 ) is necessary to use Lemma 2.6 We have 

exp z ^ 



m=l 



11/' \ '^^^ 

m=l I m=l 



. / b \ / /(A) 

EMnc...^"")^"' n( 

Ahn ^ \m=l / \m=b+l 



(n) 



(3.5) 



We now apply Lemma 

oo 



2.6 



with a„ 



n=0 



exp I 2 ^ SmC^^^ 
m=l 



^^mC*^™, if 1 < m < 6, 
6m, lib < m. 

b 



We get 



- E n 



7« sr*'") 



l(X) 
. m=b+l 



(n) 



exp 



exp 



6 e*"" 



E 

s.m=l 



m=b+l 



m 



l)r Ge(t). 



(3.6) 



This proves (3.1). The proof of (3.2) is very similar to the proof of (3.3). The only 
difference is that one has to apply Theorem 2.9 if ge{t) € J^{r,6) and Theorem 2.12 if 
ge(t) £ eJ-'(r, 6^, 7). We thus only prove (3.2). 

The function J2L=i ^(e^"" - l)t™ is a polynomial in t and is therefore holomorphic on the 
whole complex plane. We can thus apply Theorem 2^ with 

On 



S(t) = S(t,w) = exp y —(e''"^ - l)t 



k 771=1 



to obtain: 



h„E 



exp h 5^ SmC!-^ 



m=l 



Sir 



r"r(0) 



o 



n 



1-2 



(3.7) 



For ge{t) G eJ^{r,9^'~)) we get a similar expression. We have computed in Lemma 2.13 that 



1 ,n m 

f 

hr, e^n^~^ ' 



This together with (3.7) proves the theorem for ge{t) G J-'{r,6). The argumentation for 
ge(t) £ eJ-'(r, 9, 7) is similar and we omit it. □ 



4. The total number of cycles 

We prove in this section a central limit theorem and mod-Poisson convergence for -ft'on- From 
the mod-Poisson convergence we deduce Poisson approximation results for Kon as well as 
large deviations estimates. As before we use generating functions. 
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Lemma 4.1. We have for each w & C as formal power series 



n=0 



n=0 

If geit) is of class F{r,6), then 
Ee [exp(zsi^'on)] 



exp I w ^ Cm 

m=l 



^e(e"-l)gi^(e"-lj 



o 



(4.2) 



with 0{.) uniform for bounded w. If geit) is of class ej^{r,9,'y), then 



Ee [exp{isKon)] 



n 



\og{n) 



(4.3) 



9(e"-l) i<'(e"-l) ^(6*) 

r(^e-) 

with 0{.) uniform for hounded w. 

Proof. To prove the first part, one can use exactly the same argumentation as in ( 3.1[ ). We 
thus omit the details. To prove the second part, we use Theorem 2.9 for ge{t) G J^{r,9) to 
obtain 



r] [exp(e'"^7e(t))] 



o 



with 0{.) uniform for bounded w. We thus get with Lemma 2.13 

m 



nee' 



O 



n 



n 



Ee [exp(iisiron)] 

The argumentation for ge{t) G eJ^{r,6,'j) is similar. 
We can now prove a central limit theorem for K^n- 

Theorem 4.2. Let g0{t) be of class J^{r,9) or of class eJ-'{r,9,j). We then have 

i^on -eiogjn) ^ ^^^^ 

Proof. We prove this theorem by showing 

isKon 



(4.4) 

(4.5) 

□ 



E 



exp 



n 



Jsdy/login) -93^2 



(4.6) 



(4.7) 



and then applying Levy's continuity theorem. Since 0{.) in (4.2) resp. (4.3) is uniform in w 
on compact sets, we can chose w = e*'*/V^°§("). We get 



Ee 



exp 



Vlog(n) 



On 



~ n 



exp ^ log(n) 

isey/los{n) -03^2 



exp ( 6 \og{n) (ev'i°g(") — l) 

+ o(log-^/^( 



^/\og{n) 2 log(n) 



□ 
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In fact Kon converges in a stronger sense, namely in the mod-Poisson sense: 

Definition 4.3. We say that a sequence of random variables Zn converges in the strong 
mod-Poisson sense with parameters A„ if 

lim exp(A„(l - e^'))E \e''^-] = ^(s) (4.8) 

n— ^oo 

locally uniform for each s G M and "^{s) a continuous function with \E'(0) = 1. 



The mod-Poisson convergence is stronger than the normal convergence in Theorem 4.2 



since mod-Poisson convergence implies normal convergence (see [231 Proposition 2.4]. Mod- 
Gaussian convergence and mod-Poisson convergence were first introduced in ^23j and [21j . 
Details on mod-Poisson and mod-Gaussian convergence and its use in number theory, prob- 
ability theory and random matrix theory can be found in [23] and [21] . 

Theorem 4.4. Let ge{t) be of class J-'{r,6) or of class eJ-'(r, ^,7). Then the sequence K^n 
converges in the strong mod-Poisson sense with parameters K+9 log(n) with limiting function 

m 

r(ee") ■ 



Proof. This follows immediately from Lemma 4.1 □ 



Remark: 

In fact we could show mod-Poisson convergence with parameter 6'logn instead of i^r + 6'logn; 
then the limiting function would be slightly modified. 

It is natural in this context to approximate Ko^n with a Poisson distribution with parameter 
\n = K -\- ^log(n) (or simply A„ = 6'log(n)). To measure the distance between Kon and a 
Poisson distribution, we introduce some distances between measures. 

Definition 4.5. Let X and Y be integer valued random variables with distributions fi and 
V. We then define 

(1) the point metric 

dioc{X, Y) := d/oc(/i, ^) ■■= sup l/i {j} - V {j} I (4.9) 

jez 

(2) the Kolmogorov distance 

dKiX,Y) := dK{n,u) := sup |/i{(-oo,j]} - z/{(-oo,j]} | (4.10) 

iez 

Lemma 4.6. Let PK+eiog{n) be a Poisson distributed random variable with parameter K + 
6\og{n). If ge is of class J^{r,9) then 

dloc{Kon, PK+eiog(n)) ^ (4-11) 

log^nj 

dKiKon, PK+e\og{n)) < , . . (4-12) 

A/log(n) 

with ci > 0, C2 > independent of n. If ge is of class eJ^{r,6,'y) , then 

dloci^on, PK+e\og(n)) < — (4-13) 

log(n) 
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Proof. Formulas (4.11) and (4.12) can be proven with Proposition 3.1 and Corollary 3 .2 in 

m 



[2] for X = e 



(E'+eiog(n)){e'^-l) 



, ipfi{s) = 1 and e the error term in Lemma 



4.1 



We cannot prove (4.13) with Proposition 3.1 and Corollary 3.2 in [2], since the characteristic 
function of K^n does not fulfill the requested conditions. But we can modify the method in 
[2]. We have by the Fourier inversion formula 

1 



2n 



(0^(s) - 0j.(s)) ds 



(4.14) 



with (f)n{s),(j)u{s) the characteristic functions of fi and u. The characteristic function of Kon 
is given by (4.3) and the characteristic function of PK+eiog{n) is exp((i^' + 6* log(?T,))(e*'' — 1)). 
We thus get 



< 



< 



[Kon 
1 

2^ 
1 

2^ 
1 

2^ 



j]-P[Peiog(n)=j] 
giiSg(K+eiog{n))(e"- 



-{K+eiog{n))s^/2 



r(^) 



r(^) 



1 +0 



log(n) 



ds 



g-(E-+eiog(n))sV2^|_g| 



ds + 



ds + 



login) 



log(n) 



We are now in the same situation as in the proof of Proposition 2.1 in [2j. One can thus use 
exactly the same arguments to get the desired upper bound. We thus omit the details. □ 

We now wish to deduce some large deviations estimates from the mod-Poisson convergence. 
For this we use results from some work in progress [23] which establishes links between mod-* 
convergence (e.g. mod-Poisson or mod-Gaussian convergence) and precise large deviations. 
More precisely the framework is as follows: we assume we are given a sequence of random 
variables X„ such that '{)ri{,z) = K[e^^"] exists in a strip —e < Re{z) < c, with c and e positive 
numbers. We assume that there exists an infinitely divisible distribution with moment 
generating function exp(?7(z)) and an analytic function (f){z) such that locally uniformly in z 

exp {-tnr]{z)) (fn{z) (t){z), n -> oo, (4.15) 

for — £ < Re{z) < c and some sequence tn tending to infinity. We further assume that (f){z) 
does not vanish on the real part of the domain. 

Theorem 4.7 ([25J). Let {Xn) be a sequence of random variables which satisfies the assump- 
tions above. Assume further that these and the corresponding infinitely divisible distribution 
have minimal lattice N. Assume further that the rate of convergence in (4-15) is faster than 
any power ofl/tn. Let x be a real number such that tnX G N and such that there exists h with 
ri'{h) = X. Noting I{z) = hx — ri{h) the Fenchel-Legendre transform, we have the following 
asymptotic expansion 

-tnl{x)) 



Xtr 



exp I 



^j2Txtnr]" {h) 



m 



— + — + 
t t"^ 



, n — )■ oo 



(4.16) 



where {4-16) has to be understood in the sense that for every integer N, one has 



P[X„ = xt. 



exp(-tnJ(x)) 

^^TXtrji'ih) 



ax 



+ . . . + 



ajv_i 
'tE^ 



+ 



1 
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The coefficients (aj) are real, depend only on h and can be computed explicitly. 

We now apply the above result to the random variable X„ = Kon — 1 in order to have a 
random variable distributed on N. It is easy to see that in this case t]{z) = — 1, h = logx 



and I{z) = z log z — z + 1. It follows form Theorem 4.4 and Lemma 4.1 that Kon — 1 satisfies 



T(6) 

the assumptions of the above theorem with (piz) = — — r — and t^ = K + 9 log n. Using 

r(e^)r(6'e^) 

Stirling's formula we obtain the following large deviations estimates (the case K = and 
6 = 1 corresponds to the result of Hwang [19] ) . 

Theorem 4.8. Let X„ = K^n — 1- Let x E M. such that t^x G N, where tn = K + Ologn. 

Then 

P[X„ = xt^] = e-'-% [t^^^h^, + O 



k\ \T{x)T(9x) \logx 
Remark: 

In fact one could obtain an arbitrary long expansion in the above above theorem. 

5. Some examples 

In this section we consider some examples of sequences 9 and check if ge(t) is of class J-'(r, 6) 
or of class eJ-'(r, 9, 7). 

5.1. Simple sequences. 

5.1.1. The Ewens measure. The simplest possible sequence is the constant sequence 6^ = 9. 
This case is known as Ewens measure and is well studied, see for example pQ. We have 

g^(t)=9\og{j^^. (5.1) 

We thus have g@{t) G ^-"(1, and our results apply. 

5.1.2. The condition , \Sm.-e\ ^ have for t — 1 and \t\ < 1 

geit) = E -^"^ + E —t"- = ^log ( ) + E —t"- 
^-^ m ^-^ m \1 — t J ^-^ m 

m=l m=l ^ ' m=l 

=^iogf-^Vf:— +£— 

\l-t J ^ m ^ m 

^ ' m=l m=l 

.«,„,(_l_),f;^,„(i). ,5,) 

^ ' m=l 

It is clear that gQ{t) G eJ^{l,9,'~f) if we assume \9m — 9\ = 0{m~'^). It is not clear from 



(5.2) if gisif) G ^-"(1,^) or not, even when one can extend geit) holomorphically to some 



A(l,i?, 0). If a holomorphic extension is available, then one can modify Theorem 2.9 by 



replacing 0{t — 1) by o(l) in the definition of ^-"(1,^). The only difference in the result is 



that one has to replace the error term O(^) in (2.12) by o(l). We do not prove this here 
since one only has to do some minor changes in the proof of Theorem 2^ 
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5.1.3. The condition Yl^=i \^rn — 6*1 < oo. We have for t — )• 1 and |t| < 1 



ge{t) = 9 log 
= 6 log 
= eiog 



1 



1 -t 
1 

1 -t 
1 

1-t 



m=l 

oo 

m=l 

oo 



m 



6m — 6 



m=l 



m 



(t*" - 1) 



m 



m=l 



m 



' m— 1 



m=l 



m 



+ 0{t-l] 



(5.3) 



As before ge(t) ^ eJ-'(l, 6*, 7) if we assume \9m—0\ = Oim^"'). If geit) can be holomorphically 
extended to some A(l, i?, 0), then (5.3 ) shows that geif) G J^(l, 6*) (if we have an asymptotic 
expansion near 1 in A(l, i?, 0)). 

5.1.4. The sequence 6^ = e~"™ and Yl'^=i ^"'^"^ < 00 or X]m=i 1*^™ ~ "^l < 00. Both 



conditions on am ensure that |e "™ — e 
as above to see that in the first case 



(7e(t) = e-°log 

and in the second case 

(?e(t)=e-"log 



< C\am — <y\. One can use the same argumentation 



1 -t 



E 



m=l 



m 



+ E 



m,=l 



+0(1) 



O(t-l). 



(5.4) 



(5.5) 



This shows that we are in the same situation as in Section 15.1.21 and Section 15.1.31 

5.2. polylogarithm. Let 9m = with 5 7^ 0. We then have 

00 ^ 

Lin-5 :=ge{t) = Y,^t"^- 



m=l 



(5.6) 



The functions Lii+5 are known as the polylogarithm. A simple computation shows that the 
convergence radius of Lii+5 is 1. It was shown by Ford in [16] with Theorem 2.10 that the 
polylogarithm can extend holomorphically to the whole complex plane split along the axis 
M>i. The asymptotic behaviour of Lii_|_5 near 1 can be computed with the Mellin transform. 
This has been done in [TTl Section 3]. We just state here the result. 
The case 5 = is trivial since we then have Lii = — log ( jz^) . 
If 5 e {1,2,3,---}, then 

Lii+, = C(5 + l) + 0(t-l) (5.7) 

and for (5^ {0,1, 2,---} 

Lii+5 = r(-5)(l - tf + ({6 + 1) + 0{t - 1) (5.8) 

For 5 > we can apply Theorem |2.9[ but the main term in the asymptotic expansion is 0. 
If (5 < then Lii+5 is not of class F{r,6) nor of class eJ-'(r, 6^, 7). This shows that we 
cannot apply Theorem 2.9 nor Theorem 2.12 Ercolani and Ueltschi have shown in [9J 
that Kq^i converge for 5 < in distribution to a shifted Poisson distribution (without re- 
normalization) and for 5 > they have shown E [-R'on] ~ An^ . It is not yet known if Kq^ 
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converge in distribution (after re-normalization). The method of Ercolani and Ueltschi bases 
also on generating functions, but they use the saddle point method to extract the asymptotic 
behaviour of [t"] [g\. In this way they can have weaker assumptions on the sequence G, but 
do not get information on the error term. 



5.3. 9m = exp(cm^). A simple computation shows that the convergence radius r of ge{t) is 

(5.9) 



1 if ^ < 1, 

e-^ if ^ = 1, 

if ^ > 1 and c> 0, 

oo if 6* > 1 and c < 



We can not apply our method for ^ > 1 and for 6 = 1 we have ge{t) = — log(l — te'^). We 
thus restrict ourselves to 6* < 1. 

Lemma 5.1. For 6 < 1, the function ge{t) can holomorphically extended to the whole 
complex plane split along the axis M>i. 



exp(c2: ) 



Proof. As simple computation shows that 
9 < 1, we can find for each A g]0, 7r[ a C > with 

0xp(c2*) 



we can apply Theorem 



2.10 



< 2exp(|c||zn for Re{z) > \. ^ 



mce 



cxp{cz'') 



with 



< Cexp(yl|2;|). This shows that 

□ 



We now determine the asymptotic behaviour near 1. 
Case g < 
We have 



fTi 

,e(t) = Y: eMcm^- = E E ^^""'^"^ = E ^Li,_.,(t). 



oo oo J, 
C 



oo J, 
C 



(5.10) 

771=1 A;=0 771=1 A;=0 

We know the expansion of each summand near 1, but we have to justify that we plug them 
in. We first look at the case |t| < 1. We use the same argumentation as in (5.7) to see that 

oo ^ 



m=l 



-kd 



oo ^ oo ^ / m \ 

E m^-ke + ('^ ~ E m'^-ke ( ) 

m=l m=l \fe=l / 

T71=l \ 77t = l / 



with 0{.) uniform for —kO > 2. We thus can this put into (5.10) and get for |t| < 1 

1 



geit) = log 



1-t 



+ K + 0{t-l) for t ^ 1, |t| < 1 



(5.11) 



(5.12) 



One can now use the same argumentation as in [Til Section 3] with the Mellin transformation 
to see that the expansion is also valid for t in some Aq. 

Case c < 0, < g < 1 

In this case it is easy to see that ge{t) = K + 0{1 ~ t) for t — )■ 1. We thus cannot apply 
Theorem 2^ nor Theorem 2.12 This case was also considered by Ercolani and Ueltschi. 
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They have shown that Kon converge as in polylogarthm case in distribution to a shifted 
Poisson distribution. 

Case c>0,0<g<l 

The behaviour of K^n in this case is an open question. 

6. The generalized weighted measure 
In this section we introduce the generahzed weighted measure. This is defined as 
Definition 6.1. Let '■ N — > ]R>o be given for m > 1 with -Fm(O) = 1. We then define 

(6,1) 



771=1 



with hn{F) a normalization constant. 

It follows immediately from the definition that Pe [.] = Pp [.] with Fjn{k) = 6^. 

Our approach has been based so far on generating functions, and more especially on Lemma 2.6 



It is obvious that Lemma 2^ cannot be used anymore for general functions Fm, but we can 
prove a more general version of it: 

Lemma 6.2. Let Am : N — > C for m > 1 be given with Am{0) = 1 for m > 1. We then have 
as formal power series 



E 7^'" n i^n^iCm)) = n EG Urn, ^) wUh EG(A,t) := ^t' 

X ^■'^ m=l m=l ^ ^ fc=0 



with Cm = Cm(A) as in Definition 2.2 
Proof. We have 



We get 



E n i^rniCm)) = ^ 



11™=! ^m(Cm)^ 



Ahn m=l 



ci,- ,c„eN 
E^=i'»c„=n 



E n 

ci,---,c„eN m=l 
Em=i"^c„=n 



Cm! \ m 



(6.2) 



(6.3) 



1 A 



A m=l 



E n 

(c^)-=i, m=l 



Am,\Cr\ 



C ' 



JJ EG (a 



m=l 



m 



by definition of EG(F,t). 
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oo I CO 

HE 

m=l \fc=0 



Am(fc) / 1 



m 



(6.4) 

□ 



As before we can write down the generating functions for the cycle numbers and for Kon. 
Theorem 6.3. We have 



n=0 m=l 

If b eN and Si, ■ ■ ■ , s;, G M are fixed, then 



m 



n=0 



exp z 



n) 



m=l 



\{eg[f, 



7Tl=l 



\{eg[f, 



m=l 



m 



and for each w E C 



^r/i„(F)E [exp (wKon)] =1[EgIf, 



n=Q 



m=l 



m 



(6.5) 



(6.6) 



(6.7) 



to prove (6.6). We have 



Proof. The identity (6.5) follows from (6.6) by choosing Si 



Sb = 0. It is thus enough 



J2 hn{F)Ep 



n=0 



exp Z ^ SmC^ 



= 5^ ^t'^l exp U J2 s^Cl^n n FUC„ 



m=l 



m=l 



i(A) 



n)^ 
m I 



n PrriCi 



n)\ 
m ) 



, m=6+l 



We now use Lemma [6.21 with 

^m(^) 



e''-^YFra{}i), if 1 < m < 6; 



Fmi}^), 

A simple computation then shows that 



if m > 6. 



EG A 



EG\Fra,e 



for 1 < m < 6. This proves (6.6). The proof of (6.7) is similar. 



(6.8) 
(6.9) 

(6.10) 

□ 



We have thus found the generating functions in this general setting. To get the asymptotic 
behaviour as in Section |4] and Section |3j we need some analyticity assumptions. We now 
give here some examples 

6.1. exp-polynomial weights. Suppose that a polynomial P(t) = Qt\ Ylk=2 ^^t^ given 
with 6fc > and > 0. We then define F^ implicit by the equation 

EG(F„,t) =exp(P(t)). (6.11) 



It is easy to see from (6.11) that Fm(0) = 1 and Fm{c) > 0. The functions Fm thus generate 



a probability measure on Sn- We do not need an explicit expression for Fm since we are only 
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interested in the asymptotic behaviour of Cm and i^on- We have for \t\ < 1: 

^ oo , /+m\\ / / °° +m\ oo 

Y^f'K = n (^jj = H ^ b r s 

n=0 m=l \ \ / / \ \m=l / k=2 m=l 

= exp |^-^^log(l - t) + ^6,Li,(t'=) j (6.12) 
with Uk{t) ■= YZ=i B- We get that -9 log(l -t) + Yfm=2 Lifc(^') is of class ^(^, 1) and we 



can apply Theorem 2.9 It is now obvious that we get in this situation the same asymptotic 



behaviour as for the Ewens measure. 



6.2. spatial random permutations. The measure in this section comes from physics and 
arises from a model for the Bose gas and has a connection to Bose-Einstein condensation. 
We give here only the definition and a very brief idea of the model behind, but avoiding 
further details. We refer to ^ for a more comprehensive overview. 

We follow here the physicists notation and write 9m = e~"™. Let ^ : M'^ — t- M be a continuous 
function and A be a lattice in be given. Assume that for all /c G A. 

0<e^^W:= I e-«(^)e-2"^<'='^> rfx and Ve-^('=)<oo. (6.13) 



fceA 



We then define 

F^{c) = Fj^'^\c) = j e-"™ ^ e"^^'^)™ j . (6.14) 
V fceA / 

The functions Fm are well defined since there exists only finitely many keA with e-"^'^)™ > 
1. As before, we use hn = hl^'^^ for the normalization constant. 

We now describe the physical model behind this. Let D be a fundamental domain of M'^/A 
and Xi, ■ ■ ■ , x„ be n particles in D. The function ^ plays the role to penalize certain config- 
urations of the Xi. More precise, the probability of a given configuration is defined as 

1 / " n \ 

^n,dx W, dx] = exp - ^ ^{Xm - X^^m)) - J]] Cm<^m ■ (6.15) 



m=l m=l 



where dx is the normalized Lebesgue measure on D. The probability measure Pn,dx[-] now 
induces a probability measure P„[.] on Sn (by averaging over D\ It is not obvious, but it 
can be proven that P„[.] and Pp [.] are the same measures, see [U Proposition 3.1]. 

We now compute the asymptotic behaviour of Cm and K^n with respect to this measure. 
The lattice in [1] depends on n and is so chosen that the density p := is fix. We assume 
here that the lattice A is independent of n. The reason why we do this that we need in our 
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approach that the weights are independent of n. We define 



-e{fc)\m 



(6.16) 
(6.17) 



fceA m=l 



The functions g^"\t) and g^'^'^^t) are both formal power series in t since by assumption 
g-<:(fc) > and the coefficient of each t"^ is finite. 

The function g^^'^t) agrees with geit) for 6 = ie^°'"')men- '^^^ difference to the Section [s] 
and [4] is that g^°''^\t) will play role of ge{t) and not g^^'^t), i.e. the asymptotic behaviour of 
Cm and Kom depends directly on the analytic properties of g^°''^^ (see below). 
We use the same argumentation as in the previous sections to get the asymptotic behaviour 
of Cm and Kq^- We begin with the generating functions: 

Theorem 6.4. Assume that X]a;ga ^"'^'''^^ convergent. We then have the following identities 
( as formal power series ) 



G("'«)(t) := = exp (^("'^H^)) 



(6.18) 



n=0 



n=0 



G("'«)(t) fjexp (e' 



-11 



m=l 



m 



e ^ e 

fceA 



-e(fc)m 



5];/i^«)E [exp{wKon)] = exp {wg'^'''^\t)) . 



(6.19) 
(6.20) 



n=0 



Proof. We begin with (|6.18|). We have 

t 



EG Fm, 



m 



E 

c=0 



F™, c ft 



cl \ m 



E 



1 f- 



exp 



e ^ e 

fceA 



c=0 
e(k)m 



cl \ m 



E' 

fceA 



-e(fc)m 



(6.21) 



Thus 



l[EGlFm, 



m=l 



m 



m 



exp Yl 

\m=l 



e ^ e 

fceA 



't{k)m 



(6.22) 



vfceA m=l 



^g-e(fc)^)m 



expK^^7"(e-^('=)t) 



vfeeA 
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The above expressions are formal power series since g^^'^^t) is a formal power series. This 



proves (6.18). We have 



EG F^,,w- 



m 



exp w — I e 

m 



fceA 



-e{k)m 



The identities (6.19) and (6.20) thus follow immediately from Theorem 6.3 



(6.23) 



□ 



The next step is to obtain the asymptotic behaviour of Cm and -ft'on- For this we have to 
assume some analytic properties. In this context, it is natural to assume that of class 

J'{a,r). 



Lemma 6.5. Let g^"^ be of class J-'{a,r) . We define 

r := min {e''^^\ k e A} and A := # {fc G A, e"^^^ 
Then g^^^'^^t) is of class J-'{Aa,rr) . 



(6.24) 



Proof. We have by assumption that holomorphic in Aq (r,i?,0) for some (f),R. Thus 

(^(")(^g-<:(^)jf;^ is holomorphic in e'^^''^Ao(r, R, 0). This shows that g^^^'^^t) is a sum of functions 
holomorphic in rAo(r, i?, 0). We thus have to show that the sum is convergent. We have 
Ig^^^t)] < Ci\t\ for \t\ < r/2 since by assumption ^'^"•'(O) = 0. Let t in rAo{r,R,(f)) be fixed. 
Since the X]a:ga convergent, there are only finitely many k G with |t|e^^'^'^'' > |r. 

We thus get 



|^7("'f)(t)| 



J2g^-\e-<'h) 

fceA 



<C2+J2 Cile-'^^'hl <C2 + Ci\t\ J2 \(^~"^^^\ < oo- (6-25) 



fceA 



This proves that g^°''^\t) is holomorphic in rAo(r, R, 0). The other statements are clear. □ 



We now apply Theorem 2.9 to obtain: 



Theorem 6.6. Let g^"^ be of class J-'{a,r) and r,A be as in Lemma 6.5 We then have for 
each 6 G N 



with Yi, - ■ ■ ,Yh independent Poisson distributed random variables with 



(rr)' 



m 



e " ^ e 

fceA 



-t{k)m 



(6.26) 



(6.27) 



Proof. The proof of this theorem is the same as the proof of Theorem |3.1| One only has to 

□ 



use Theorem 16.41 and Lemma f6. 51 We therefore omit the details. 



Theorem 6.7. Let g^""^ be of class J-'{a, r) and r, A be as in Lemma 6.5 Then converges 
in the strong mod-Poisson sense with parameters Aa\og{n) with limiting function ^(^tj^- 
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